The violation of the weak energy condition, Is it generic of spontaneous scalarization ? 
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It was recently shown by Whinnett & Torres Q that the phenomenon of spontaneous scalarization 
(SC) in compact objects (polytropes) was accompanied also by a spontaneous violation of the weak 
energy condition (WEC). Notably, by the encounter of negative-energy densities at several star 
places as measured by a static observer. Here we argue that the negativeness of such densities is not 
generic of scalar tensor theories of gravity (STT). We support this conclusion by numerical results 
within a class of STT and by using three realistic models of dense matter. However, we show that 
the "angular parts" of the additional conditions for the WEC to hold pcff + T^^°^^ > Q tend to be 
"slightly violated" at the outskirts of the star. 

PACS numbers: 04.50.+h, 04.25.Dm, 97.60.Jd 



^ ■ I. INTRODUCTION 

> 

ps| \ Scalar-tensor theories of gravity (STT) are one of the simplest alternatives to Einstein's general relativity (GR) [31 ■ 
T— I A particular example of these, is the Brans-Dicke theory. The general feature of STT is the introduction of a new 
fundamental scalar field which couples to gravity in a non-minimal fashion (NMC). Nonetheless, the NMC is such 
that it preserves the Equivalence Principle, because the scalar-field couples to the ordinary matter only through the 
spacetime metric. 

The NMC in STT unlike simple Einstein-Higgs systems, gives rise to field equations that can differ substantially 
^ ' from the usual Einstein's equations. Therefore, when applied to the astrophysical and cosmological settings these 
y ' theories can produce severe deviations from the general relativity predictions. The fact that STT introduce in general 
5-H . only a few new parameters make them in principle easily testable. Perhaps the most notable probes in this direction 
are those connected with the binary pulsar, the primordial nucleosynthesis, the Cosmic Microwave Background, and 
^ the luminosity distance with Type-II Supernova. The success of GR lies in the fact that with suitable matter models, 
the theory is compatible with the available information in the above systems. Modifications of the gravitational sector 
in the way of STT lead thus to a serious change in the geometry-matter relation that can upset the compatibility 
' of the theory with the observations. The stringent bounds on STT make difficult to accommodate these theories to 
all the observations while at the same time predicting new observable phenomena. One of the potentially observable 
new phenomenon is precisely the phenomenon of spontaneous scalarization (SC) arising in compact objects. As it 
was first shown by Damour & Esposito-Farese STT can induce non-perturbative effects in neutron stars (NS) 
which, roughly speaking, consists in a phase-transition that endows a compact star with a new global quantity, the 
scalar charge. This "charge" is the analogous of the magnetization in ferromagnets at low temperatures and the 
central-energy density (or equivalently the baryon mass) of the object plays in turn the role of the temperature in 
spontaneous magnetization. Damour & Esposito-Farese also shown that the consequences of SC can be dramatic 
even if the boundary-conditions (asymptotic conditions) for the scalar-field are chosen so as to satisfy the solar-system 
bounds. In particular, they showed that the binary-pulsar dynamics is very sensitive to the new phenomenon, and 
therefore, that STT can be further constrained even if they successfully pass the local and the cosmological tests. 

In this paper we will not be concerned with the confrontation of STT with the experiments, but rather with a 
result obtained by Whinnett & Torres who presented evidence that the the phenomenon of SC in NS is linked 
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to a violation of the weak energy condition (WEC). Notably, Whinnett & Torres Q found that the effective energy 
density of STT as measured by a static observer and defined as Peff n°-n^Gab/ {S>t^) (where Gab is the Einstein 
tensor and n°- a time- like vector parallel to the static Killing vector is negative and suggested that such a finding 
might be independent of the class of STT used (although depending of the chosen value of the NMC constant of the 
theory). A negative peS is a sufficient condition for the violation of the WEC. The authors argue that the violation 
of the WEC in NS might led to their instability (even for masses quite below to their corresponding Landau- Tolman- 
Volkoff-Oppenheimer limit), a very undesirable feature that would presumably rule out the STT as an alternative for 
a viable spacetime theory. 

While it is true that the effective energy density Peff of STT as defined above is not automatically positive definite 
due to the presence of second order derivatives of the scalar field, the main goal of this paper is to show that there 
are classes of STT that do not violate the condition p^g > for the problem at hand. As we shall argue in Sec. Ill, 
the violation of the condition pcff > as found by Whinnett & Torres is linked to the use of a parametrization d 
la Brans-Dicke which prejudices the choice of the specific theory which can potentially violate such condition. In our 
case, we use a parametrization of the STT where the effective gravitational constant can be chosen to be positive- 
definite. In such a case, as it is shown in Sec. Ill, the theories parametrized in both ways may be very different (the 
transformation that would be necessary to change from one parametrizaction to the other might not be well defined) . 

Altough the class of STT analyzed by us satisfy the condition pcs > for the NS of interest, they are prone to 
violate "slightly" the angular parts of the additional conditions Pcff + T^j'^^'^ > for the WEC to hold. Since this 
violation is rather weak and occurs when the energy-density and pressure of the star fluid is very low (near the star 
surface), we argue that the stability of the neutron star is not in jeopardy, and therefore that a large class of STT 
cannot be considered as unviable spacetime theories only by this fact. 

II. THE MODEL 

The general action for a STT with a single scalar field is given by 

S = J |Y^i^(<^)i?- Q(V0)2+l/(0)^|v^d4x + 5„,att . (1) 

where the matter sector will be represented by a perfect fluid. 
The field equations obtained from the above action are 

Gab = 87rT°f , (2) 
T'a! = G,s{T[, + T^, + Tr'') , (3) 
T[b = ^[Vaid^FVb^)-9abVcid^FV'<f>)] , (4) 



T^at = (Va0)(Vfc0)-g„b 



(5) 



TT'' - (P + p)UaUb + CJabP , (6) 



□ = 



Fd^V - 2id^F)V + ^{d^F) 




[1+ ) 




F+' 





(7) 



where T^att stands for the trace of T,^,ttt' Gab = Rab - ^abR, and G^s = y. 

It is to be emphasized that the action corresponds to STT in the physical (Jordan) frame, as opposed to the 
one written in terms of an unphysical (conformally transformed) metric (the Einstein frame). Hereafter we consider 
F(0) ^0. 

In a previous work j^j (hereafter SSN), we analyzed the SC phenomenon in detail for the class of STT given by 
F(0) = 1-1- 167r^(/)^ (here ^ stands for the NMC constant) and for spherically symmetric (static) neutron stars. In that 
analysis, we exhibited the SC phenomenon for the cases ^ = 2,6 and employed three representative realistic equations 
of state (EOS) for the nuclear matter: the model of Pandharipande 6] (hereafter PandN; representative of a "soft" 
EOS), the model II of Dfaz-Alonso Q (hereafter DiazII; representative of a "medium" EOS), and the model "0.17" 
of Haensel et al. (hereafter HKP; representative of a "stiff" EOS). Our conclusion in SSN was that the SC is a 
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phenomenon that ensues independently of the details of the EOS, but that strongly depends on the compactness of 
the star (high central energy-density). Moreover, we provided a heuristic analysis that clarifies why NS configurations 
with SC are energetically preferred over those with a trivial scalar field, despite the unintuitive expectation coming 
from the fact that the effective gravitational constant GcS = 1 /F decreases as the scalar field departs from the trivial 
configuration, which in turn leads to a reduction of the absolute value of the star's negative binding energy. 

In order to analyze the SC phenomenon we computed hundreds of NS configurations. Now, to inquire about the 
status of the WEC, we take representative data of our catalogue and recompute the energy-density profiles Pcsif) 

as well as the other two relevant quantities Pcs + T\^''''^ and Peff + T^g (where the spherical symmetry implies 



T 



0(off) 



T^^""^^). The WEC will be violated if any of the three above quantities is negative. 



For the static and spherically symmetric metric gab — diag(— A^^, A^, r^, sin^ 0), the effective energy-density 



Peff := n'^n^T^f and the quantities peff + Tr and pcS + T g ' , respectively [where T^f is given by Eq. and 
stands for the timelikc unit vector associated with a static observer] we have 



Pc 



PcS = 



I rjir (efF) 
B + J- r 



PcS + T' 



e(off) 



Goff 



1 + 192^e'0'Gcff 
+p + 64^e'0'Geff (2p + 3p)] , 
Gcff 



1 + 19277^2026-^(1 
1 



1 drN 



\r N 



(1 + 192^f 02Gcff) 



+ j^{dr(^f (1 + 4e -f 1287rC202Gcff) +p + p+ 128t:^^ cj)^G,s (p + 3p) 



G. 



cS 



A^(f>dr(f) I drN 
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+ P + P 



(8) 

(9) 
(10) 



which explicitly shows the usual contributions of the the perfect- fiuid (p and p) plus those arising due to the NMC. 
Here N — N/Nq is the renormalized lapse with respect to its value at r — 0. Clearly when (f>(r) = , one simply 

recovers poff = P, and peS + ^ r = p + P = PeS + g ■ 

We have integrated the system of equations with adequate regularity and boundary conditions (leading to asymp- 
totically flat spacetimes) and with 0(r oo) ~ Qs/r where Qs is the scalar "charge." For a fixed ^, all the NS 
models are parametrized by the central energy-density of the fluid pq. Therefore a critical value Pq^^^ (or equivalently 
-^bar ) characterizes the onset of SC. Figures depict pcs (left panels) and PcS + T^r ''^^ (right panels; sohd curves) 
and pcs + T^^°^'^ (right panels; dashed curves). The lower mass {M"^^) curves of Figs. indicate the critical 

configurations towards SC while the rest of the curves correspond to NS endowed with a "scalar charge" Qs- 

The important point to note is that none of the curves of Figs. ^-J^show a violation of the condition p^s > 0. From 



the figures we appreciate that the condition p^s + T r > is also satisfied. However, the quantity p^s 
start becoming "slightly" negative near the surface of the star. 

As we have stressed above, pcS bas not an a priori definite sign, and in fact, the term linear in drcj) could be 
negative within the NS. Moreover, one can define an energy-density that is useful to understand the appearance of 
the phenomenon of SC, and is the one given as follows 0, 



e (off) 



Gcff 



1 + 1927re2</.2Gcff 



[p + M^e4>''G,s (2p + ip)] 



(11) 



which heuristically encompasses the contribution of the NMC to the energy-density of the perfect fluid. Note, however, 
that PcS = p + p^ + ■■■ with p + p^ > 0. Therefore Pcs > 0, since typically | — (j){dr(f>){drN)/ {NA"^)] ^ p within the 
NS. In our heuristic analysis (SSN) which was also confirmed by the numerical solutions, we showed that the negative 
contribution p^ more than compensates the decrease in the absolute value of the binding energy due to the reduction 
of Gcff , resulting in an overall lower total energy [the ADM-energy that is obtained from the integration of Eq. ((SJ] as 
compared with the energy of the configuration with the same baryon-mass but with a null scalar field (a configuration 
in pure GR). 
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III. DISCUSSION 



In order to give an insight for the non-negativeness of the effective energy-density let us consider eq. Q without 
specifying and compute Pcs = n'^n'^T^^: 



Pes - J 



Stt 



(V0)^ , {d^Ff dl^F 



-p+%^(2p + 3p) 



(12) 



This reduces to Eq. JHJ for F(0) = 1 + 167r^(/)^ and for the static and spherically symmetric case. 

For usual nuclear matter (p > and p > 0) and for F > 0, d'^^F > 0, only the first term is not positive semi-definite. 
Whinnett & Torres {Tl have found a negative poff at various star positions, notably at the center of the star (r = 0). 
So let us first focus on poff at r = 0. By regularity at the origin, dr4>\r=o — 0, therefore 



Poff — TT 



{d<t,F)a In 

i6^F„ y^po 



^(94'F)l 
1&-kFo 



(13) 



We can now investigate in which situations p^^ can be negative pTl |. 

a) STT with F{(f)) > (in particular Fq > 0). Clearly for this case, p'^g > (since we consider only po > and 
Po > 0). This includes the class of STT we analyzed numerically F{(f)) = 1 + WirS^cf)^ , ^ > 0, for which we did not 
encounter negative effective central energy-densities [cf. Figs. ^-13 (left-panels)]. 

b) STT with F{(j>) < for some (p (in particular Fq < 0). Clearly Fq < is a necessary condition for < 0. One 



1671- _Fo 



(2po + 3po) > and 1 



167rFo 



> 0. 



can then obtain pj^g < if in addition the following inequalities hold po 
From this, we have the following two possibilities: 1) If po — 3po > then the two inequalities hold provided 
167r|Fo| > 3{d^F)l; 2) If Bpo-Po > then the two inequalities hold provided 167r|Fo| > C{d^F)l where C = 2+3po/po- 

According to 'b)', it is in principle possible to obtain p'^g < for Fq < 0. However, this case is completely 
pathological: on one hand at the center of the star Fq is negative (which in turn implies a negative effective gravitational 
constant = 1/Fq <0). On the other hand, asymptotically F — G^^ ^ 1 > (i.e. Gcff approaches the Newtonian 
value) in order to pass successfully the solar-system bounds. Therefore F{r) — F(0(r)) must interpolate between the 
negative and the positive value, implying that at some spherical shell (inside or outside the star) F — and thus 
Gcff — > cxD, where the equations become singular. 

In order to avoid this kind of pathologies within a neutron star model one should consider only the cases F{(j>) > 
which lead us to the conclusion that any physically viable STT will not produce a negative central effective-energy 
density 

Now, although we have just proven rigorously that Fq > implies p°ff > 0, we still have to proof the same for Peff (r) 
given by Eq. ifT^ . Although this could be a difficult task, the empirical numerical evidence shows that if p^g > and 
for a class of STT (e.g. a class of STT with F > and d^^F > 0) then Pes{r) > for the same class due to the fact 

that \{d^F)n''n''VaVb(l}\ « p within a NS. 

The above heuristic analysis which is supported by our numerical results leads us to the following conjecture: 
Conjecture: Consider a STT with F(<p) > and d^^F > (in particular Fq > at the center of a neutron star 

which implies p'^g > 0), such that asymptotically F ^ 1 and with V{<p) = 0, then Pcsir) > for all the static NS 
configurations, including all those endowed with Qs (where the eguality holds only outside the star and for a trivial 
scalar field) 

Now, in order to give some insight about the violation of the WEC found by Whinnett & Torres Q , let us consider 
their Eq. (5) (where peS is denoted by p) and focus on the central energy density (where the regularity conditions 
and the spherical symmetry implies that the gradient contributions become null): 



Po 



(3po - Pq) 



2t^o($o)-H3 



(14) 



This expression is in fact equivalent to our Eq. H13|) , which we showed to be positive definite for Fq > and for usual 
neutron matter (po > and pq > 0). The translation from the Brans-Dicke parametrization used in Eq. H14() to ours 
is provided by 



(15) 
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On the other hand /io might not look as having a definite sign. We have exphcitly 

^0 = 7^- [Po(l - c^o) + SfToPo] • (16) 

where a{^) = 1/(3 + 2uj). Then, clearly due to the Brans-Dicke parametrization one could consider situations where 
CTo > 1 or (To < 0, which depending on the values of po and pq can lead to /io < 0. However, we note that the condition 
F{<j)) > (in particular Fo > 0) leading to > (or equivalently /io > 0), implies $ > and also < ct < 1/3: In 

fact, from a{^) = 1/(3 + 2uj) and Eqs. ^ one obtains, cr($) = ^(^gjpyfjie^F (thus F(0) > implies < ct < 1/3). 

Moreover, since $ — F{<f>) (assuming F{(l)) ^ const.) the inverse mapping is = / ^J^^^^^d^- Therefore the 

conditions $ > and < ct < 1/3 are required for this mapping <i> ^ </) to be well defined (the separate case ct = 
corresponds to i^(0) = const, from which one simply recovers GR). 

Hence, the restrictions on ct bounds in turn the choice of the Brans-Dicke function u>{^^). Any choice of a;($) leading 
to a violation of the condition < ct < 1/3 with <I> > will imply a lost of the connection with the original theory. It 
is in this sense that the use of the Brans-Dicke parametrization prejudices the kind of STT (those STT which allows 
values for ct to be outside the bounds < ct < 1/3) that can lead to violations of the positiveness of p^s- 

One could then use the SC phenomenon in NS to classify the STT in two types: Type-I) STT with ct($) such that 
< CT < 1/3 and <& > (which in particular arc valid at the star's center and therefore that coincide with the case 
"a)" analyzed above). Type-II) STT with ct($) (in particular ctq) not necessarily limited in the range [0,1/3) and 
with $0 > 0. The first kind of STT can be put in correspondence with the original parametrization {F{(f)) > 0) and 
these theories lead to < p^jf = /xq > 0, regardless of the EOS for nuclear matter (assuming p and p non-negative). 
The type-II is a STT which cannot be put in correspondence with the original parametrization (in terms of (j)) since 
the transformation as shown above is not well defined. The type-II theories can lead to situations where pq < 
but are not even related to the theories which in the original parametrization give rise to < (case "b)" above) 
which require Fq < 0. For instance, it is possible to have pq < with $o > with a suitable cto, as in the following 
situations: A) po/(Po ~ 3po) < ctq with 3po < Po! B) ctq < with po/(3po — Po) < Wa\ and po < 3po- Whinnett 
& Torres ,1] considered two classes of STT. The first one with ct — 2Kln<I' and k > 3. This example can belong to 
Type-IIA (and produce po < 0) if for instance ctq ~ 2, po < Po/6 and 1 < $o < 1.4. These conditions could be met 
in neutron stars with a low baryon mass like some of the ones computed by Whinnett & Torres 0. 

The second class of STT used by those authors is one with ct = — 1) and A > 5.35. This class is of Type-IIA 
or Type-IIB for <i>o > 1 or < $0 < 1, respectively. As mentioned above, the condition 3po < po in Type-IIA is met 
usually in NS with a low baryon mass. Since for this second class of STT, po < was found by Whinnett & Torres 
in low baryon mass configurations, presumably $0 > 1 there. 

We have thus given a heuristic explanation for a negative po in both of their classes of STT. 

Of course, STT of Type-II can lead to positive po's without implying that the WEC will not be violated somewhere 
within the star (see figures of Whinnett & Torres 1]). On the other hand, the fact that STT's of Type-I lead to 
positive definite po does not imply that the WEC is not violated either. However, our numerical analysis shows 
that this type of STT (theories with po > 0) indeed satisfy the condition p{r) > (where the equality holds only 
outside the star and for a trivial scalar field) and this is the motivation for our conjecture above. Along the same 
line of reasoning one could further conjecture (conjecture 2) that given a fixed EOS, if a STT of Type-II violates the 
condition p{r) > somewhere within a star (not necessarily at the center) then the STT will violate the WEC at the 
center for some other po (defined as pg ). That is, the second conjecture states that po = f{po) has always a negative 
branch (/(po^) < 0) in the Type-II theories. The first conjecture on the other hand roughly states that a STT (with 
the given assumptions) that does not violate /io > (Type-I with /io(po) > 0) will not violate it elsewhere within the 
NS (i.e., p{r) > within the star). 

We can now turn the attention to the condition poff + T^J which the numerical analysis shows to be negative near 

the star surface. From eq. pU|) we appreciate that beyond the compact support of the nuclear fluid, poff + T^J ^ 

Gcff ^^^^ - ^) . If we roughly approximate A^^ ^ N ^ {1 - 2M/r)i/2 for r > i? (where R and M stand for 

the star radius and total mass) then, the factor (^^ — ~ ^/r [(1 — 3Af/r)/(l — 2M/r)] which is positive even for 

the maximum mass configurations. Therefore p^s + T^^°^^ ^ q ^^ Ai,4>dr<P{i-iM / r) ^ gj^^^g ^j^g q^q configurations have 
dr<j) < near and beyond R [(l){r) ~ Qs/t asymptotically, with Qg > 0], then pcff + T^^'^^'' is negative (for > 0) 
near and beyond the star surface. Note from Eq. ||5J), that the term p^s + T^r ''^^ [unhke pcS + T^g '^^^; cf. Eq. l(TU|l ] 
does not become negative because of the sign difference in the term with ^(j)dr(t>, and the remaining terms involving 
the scalar field are positive definite (for ^ > 0). 

In summary, there exists classes of STT that give rise to the phenomenon of spontaneous scalarization with and 
without a violation of the condition p^ff > 0. Examples of STT like the ones we analyzed numerically does not violate 
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that condition. However they do violate the WEC "sUghtly" because peS + T g ' starts becoming "shghtly" negative 
from the star surface. On the other hand, there are STT Hke the ones analyzed by Whinnett & Torres [l| which do 
violate the condition p^s > which straightforwardly leads to a violation of the WEC. In our case it is unclear whether 
the small violation of the WEC by the STT considered here will necessarily lead to the instability of the NS. In Ref. 
@ , we have shown that the energetically preferred NS configurations (with a fixed total baryon mass) are those with 
lower ADM mass which correspond to scalarized NS as opposed to the corresponding NS in pure GR. The curves of 
the ADM as a function of the central energy-density po suggest that configurations with ADM mass smaller than the 
maximum mass models (corresponding to p™^^) are stable with respect to radial linear perturbations j5j]. Moreover, 
numerical non-linear analysis strongly suggest that ordinary NS (with weak or zero scalar charge) can decay to their 
corresponding more stable strong scalarized configurations |9| and furthermore that among these scalarized NS only 
the configurations with po > p™^'^ are unstable and can collapse into a Schwarszschild black hole like ordinary NS, by 
radiating away all the scalar-field 10]. 

If found that only the scalarized NS which violate the WEC through the violation of Peff > within the star are 
unstable, then the STT allowing such violations would be ruled out. On the other hand, it is well possible that small 
violations of the WEC through the term poff + T^g"^-* in the outskirts of the NS while respecting peff > within the 
star might not lead to instabilities and therefore that the employed STT is a viable candidate for a spacetime theory. 



IV. CONCLUSION 



Our main conclusion is that the large violations of the WEC deep inside the neutron stars as an additional feature 
to the phenomenon of spontaneous scalarization is not generic of all the classes of scalar tensor theories of gravity as 
the work of Torres & Whinnett work might be thought to suggest. We have argued that there are classes of STT (at 
least a subset with F{(f>) > and d'^^F > 0) where pcS > and where the violations of the WEC in neutron stars 

occur only near and beyond the surface of the star via the angular parts of pcS + T^^'^^^ . Strong numerical evidence 
supports this conclusion. 
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FIG. 1: 

Neutron star models constructed with the EOS PandN within the scalar-tensor theory F{(j>) = 1 + 167r^<j!>^, with ^ = 6. Left 
panel: effective energy-densities pcs- Right panel: p^s + T^r"^^^ (solid lines) and p^s + T^^"^^^ (dashed lines). In both panels 
the curves are labeled by the corresponding baryon-mass (in solar mass units). The lowest baryon-mass configuration mark 
the onset of spontaneous-scalarization. The figure also shows the maximum mass models. Here pnuc = 1.66 x lO^^kgm"''. 
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FIG. 2: 

Same as fig. [Hfor the EOS DiazII. 
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FIG. 3: 

Same as fig. [T]for the EOS HKP. 



